Hole depletion and localization due to disorder in insulating PrBa2Cu307_5: 

Compton scattering study 
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The (mostly) insulating behaviour of PrBa2Cu307_i is 
still unexplained and even more interesting since the occa- 
sional appearance of superconductivity in this material. Since 
YBa2Cu307_i is nominally iso-structural and always super- 
conducting, we have measured the electron momentum den- 
sity in these materials. We find that they differ in a striking 
way, the wavefunction coherence length in PrBa2Cu307_i be- 
ing strongly suppressed. We conclude that Pr on Ba-site sub- 
stitution disorder is responsible for the metal-insulator transi- 
tion. Preliminary efforts at growth with a method to prevent 
disorder yield 90K superconducting PrBa2Cu307_a crystal- 
lites. 

PACS numbers: 74.72Jt, 74.72Bk, 74.62Dh, 78.70Ck 



The peculiar case of PrBa2Cu307_5 stands out in the 
RBa2Cu307_5 family (R for rare earth or Y) because 
it is the only easily fabricated member which until re- 
cently was consistently insulating in the normal state and 
non-superconducting (for a review see |],|). Supercon- 
ducting samples have since been fabricated ||^,^ , some of 
which show metallic behaviour in the normal state, akin 
to that seen in YBa2Cu3 07_5. The non-superconducting 
behaviour has been attributed to various causes: (1) +4 
valence for Pr leading to hole filling in the Cu-0 lay- 
ers, (2) hybridization between the Pr 4f and O 2p or- 
bitals 0, leading to hole locaHzation in these layers, (3) 
a chain-based model for superconductivity and magnetic 
pair-breaking induced by Pr on the Ba site I). The first 
possibility has been invalidated on the grounds that Pr 
on the rare-earth site is essentially trivalent [^jsj while 
the last is unlikely if only because it fails to explain why 
non-superconducting PrBa2Cu307_5 should be insulat- 
ing in the normal state. 

The Fehrenbacher-Rice model of hole-localization in its 
various forms |^,^ remains a possible explanation, and to 
get further insight into this matter we have measured the 
electron momentum density in both YBa2Cu307_5 and 
PrBa2Cu307_5 using Compton scattering of high energy 
X-rays. This is the first high statistics and high reso- 
lution measurement in these materials, which are diffi- 
cult candidates for such experiments given the problems 
related to fluorescence background from heavy atoms 
and low count-rates due to absorption. These problems 
have been overcome thanks to high energy synchrotron 
radiation and dedicated instrumentation at 1D15B at 



the ESRF 10|. Our measurements show that the ex- 
perimental momentum density is similar to that given 
by LDA band-structure calculations for YBa2Cu307_5. 
However, in PrBa2Cu3 07_5, a striking difference is seen 
between measurement and calculation. Simple models 
show that Pr on Ba site defects (labelled Prsa hence- 
forward) and a lack of hole doping in the Cu-0 planes 
probably lead to this difference. These facts are inter- 
preted in the light of other experiments ]7|,p|,pj|-p^ 28 - 3^] 
in the Yi_xPi'xBa2Cu307_5 system. 

The experiments were performed on crystals grown 
by the flux method in the non reactive crucible ma- 
terial BaZr03, known to give uncontaminated, high- 
purity samples The YBa2Cu307_5 crystals were 
annealed under 1 bar of oxygen at 530° C for 170 hours 
yielding a Tc of 91.3 K and an oxygen content of 6.89 
0. The PrBa2Cu307_5 crystals were annealed under 
1 bar of oxygen at 500°C for 200 hours yielding an oxy- 
gen content of 6.92 [l5|] , and were insulating and non- 
superconducting. Both exhibited uniform twins in the 
a-b plane, as is usually observed after such a treatment. 
These twinned crystals were used for the experiment. 

Compton scattering, or inelastic scattering with very 
high momentum and energy transfer is a probe of 
the ground state, one-electron properties of the system 
p^ , p^ . The measured quantity J{pz), the Compton pro- 
file, is a projection onto the z-direction (parallel to the 
scattering vector) of the three dimensional electronic mo- 
mentum density p(p): J{pz) = p{p)dp^dpy. Comp- 
ton profiles were measured using the high-resolution 
scanning spectrometer at ID15B. A cylindrically bent fo- 
cusing Si(311) Bragg monochromator provided an inci- 
dent energy of 55.91 keV for YBa2Cu307_5 and 57.78 
keV for PrBa2Cu3 07_5 respectively. The radiation scat- 
tered at an angle of 159 degrees was analyzed by a cylin- 
drically bent Ge(440) analyzer. For both compounds 
three twinned single crystals were stacked to use the full 
beam size of 0.2x5 mm. Compton profiles were measured 
for the [100/010] and [110] directions, in symmetric re- 
flection geometry at room temperature with high statis- 
tics (5.10^ counts in a 0.05 a.u. bin at Pz = 0; 1 atomic 
unit (a.u.) of momentum = 1.89 A^^ ). The effective res- 
olution was 0.15 a.u. (full width at half maximum). The 
incident beam was monitored using a Si PIN diode and 
this was used to normalize the scanned proflles. The raw 
profiles were further corrected for geometrical effects, for 
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absorption in air and in the sample and for the analyzer 
reflectivity. The background was measured separately 
and subtracted. Finally the profiles were symmetrized 
and normalized appropriately. 

In the following we shall discuss the anisotropy of the 
Compton profiles, that is, the difference between Comp- 
ton profiles measured in two crystallographic directions 
(J [100/010] ~ J[iio])- This procedure has the advantage of 
getting rid of the isotropic contribution of the core elec- 
trons as well as residual background. Since both mea- 
sured directions are in the a-b plane, the structure in 
the anisotropy originates essentially in this plane, but in- 
cludes contributions from all entities of the unit cell, the 
Cu-0 chains, the Ba-0 and Y/Pr planes as well as the 
Cu-0 layers. However, from earlier momentum density 
measurements as well as other experiments Jl7|-[l9t it 
has been established that the Cu-0 chains are identical 
in both materials, being locally metallic with a corre- 
sponding Fermi surface (even in PrBa2Cu307_5). Thus, 
any differences between the two originate in the remain- 
ing entities, the Cu-0 layers or the Ba-0 and Y or Pr 
planes. 

For comparison with experiment we have computed 
theoretical Compton profiles based on LDA band- 
structure. A FLAPW 1^ calculation was performed 
for YBa2Cu307_5. For PrBa2Cu307_5, earher studies 
have shown that two of the three Pr-f electrons can ac- 
tually be considered as core states since treating 
them as band electrons does not modify the density of 
states significantly ||2^ , implying that theory predicts lit- 
tle change with respect to YBa2Cu307_5. The expected 
difference in the anisotropy in going from YBa2Cu307_5 
to PrBa2Cu307_5 is due to the hybridization of Pr-f elec- 
trons with 0-2p electrons which even in LDA calculations 
produces some hole depletion in the Cu-0 layers simi- 
lar to the Fehrenbacher-Rice model However, LDA 
calculations are known to work reasonably well only for 
the metallic phase and not for undoped insulating Cu- 
O layers in these materials. To estimate this difference 
due to hybridization we used the LMTO method to 
calculate the anisotropy for PrBa2Cu307_5 (with the Pr- 
f electrons considered as valence electrons), from which 
the LMTO anisotropy for YBa2Cu307_5 was subtracted. 
This difference should give us an upper bound on the ex- 
pected change since LDA is known to overestimate the ef- 
fect of hybridization of localized f-orbitals (see for exam- 
ple (ID). This difference was then added to the FLAPW 
anisotropy for YBa2Cu307_5 to give an estimation of 
the theoretical anisotropy for PrBa2Cu307_5. Liechten- 
stein and Mazin |^ have included Coulomb correlations 
in their LDA+U, LMTO calculations for RBa2Cu307_5 
and find that in R=Pr, with respect to R=Y, an addi- 
tional band crosses the Fermi level trapping holes and 
changing their nature, but with negligeable impact on 
the Compton profile according to our estimations. 

The anisotropy is shown as a percentage of the peak 
value of the Compton profile. The amplitude of the 
anisotropy in the theory is about 40% larger than in the 



experiment. This tendency is commonly observed, even 
in simple metals, and attributed to correlation effects 
under-estimated by LDA ||2^,^. However the shape of 
the anisotropy, which is what we use in our analysis, is 
generally very well described by LDA, even in highly cor- 
related electron systems containing f electrons, such as 
CeCu2Si2 ll^. In the figures shown we have scaled the 
theoretical anisotropy to the experimental one by a fac- 
tor of 1.4 to ease comparison. Fig. |l| shows the measured 
and calculated (solid fine) anisotropy for YBa2Cu307_5. 
The two are remarkably alike, showing that the overall 
description of the electronic structure is satisfactory. We 
also show the calculated (as detailed above) anisotropy 
for PrBa2Cu307_5 (dashed line), which, despite Pr-0 
hybridization, shows little change with respect to the 
YBa2Cu307_5 calculation. Fig. ^ shows the measured 
anisotropy for PrBa2Cu307_5 which is strikingly differ- 
ent than the one predicted by theory. 

When PrBa2Cu307_5 is synthesized it is known 
that Ba-rich or Pr-rich, off-stoichiometric phases 
(Pri-|-xBa2^xCu307_5) may be produced. There is how- 
ever substantial evidence that even in globally stoichio- 
metric samples substitution disorder concerning the Pr 
and Ba sites (spinodal decomposition for example) exists 
locally P,p8|-p0[ , especially if growth is followed by slow- 
cooling. Such disorder, due to the large size of Pr ions 
has also been observed in the case of similar sized Nd ]3C| ] 
and La 1^,0 ions. The amount of this substitution has 
been estimated to be a maximum of about 8% ||^,|2^. If 
such disorder were present, one effect would be to dimin- 
ish the contribution of the Pr and Ba-0 planes to the 
total anisotropy due to the absence of long-range order 
in these planes. Since accounting for such disorder in the 
theory is not trivial, for qualitative analysis we have used 
a simple model within the LMTO method, where we have 
neglected the contribution of Pr and Ba-0 planes to the 
total anisotropy (solid line in Fig. ||). In the same fig- 
ure we also show a molecular orbital calculation (dashed 
line) of the anisotropy due to a hole in the Cu3dx2_y2 
state corresponding to Cu-0 layers in undoped insulat- 
ing RBa2Cu3 06. Again we are limited to a qualita- 
tive approach because LDA predicts an incorrect metallic 
ground-state for RBa2Cu306 with anisotropy similar to 
that for RBa2Cu307_5. The calculation accounting for 
disorder mimics the strong suppression of the first peak in 
the measured anisotropy for PrBa2Cu307_5 which also 
tends towards the behaviour calculated for undoped Cu- 
O layers. We argue below how Prsa disorder could in fact 
lead to a situation where carriers are not doped in to the 
Cu-0 layers, effectively killing both metallic behaviour 
and superconductivity in PrBa2Cu307_d-. 

The charge-transfer model for the effect of oxygen dop- 
ing in YBa2Cu307_5 |32 3^ is well established, and de- 
scribes the mechanism of transfer of positive charge from 
the Cu-0 chains to the Cu-0 layers as oxygen is added to 
the chains. Two possible mechanisms may exist to block 
or modify this charge transfer. The first that we propose, 
supposes that if Pr is present at the Ba site, the needed 



2 



charge is supphed by Prsa (possibly in a 4+ oxidation 
state because the chemical environment on this site is 
different), inhibiting hole creation in the layers. How- 
ever, the chains would maintain their microscopic metal- 
lic character as has been experimentally observed and 
mentioned above. It is significant that according to Cava 
et al. 1^ the total charge transfer associated with Tc go- 
ing from OK to 90K is 0.08e/Cu-O layer in the unit cell. 
Knowing that the concentration of Prsa defects is a few 
percent, we have here a hole-depleting mechanism which 
is of the same magnitude as the hole-creating mechanism 
of charge transfer. Experimental observations of hole- 
depletion in the Cu-0 layers include optical conductivity 
@, NMR g. Hall-effect |3|] and transport |3^] mea- 
surements. Most experimental evidence favours a valence 
for Pr close to +3 (compatible, in our model, with Pr on 
the rare-earth site) but also indicates the probable pres- 
ence of a small fraction of -1-4 valent Pr H,^,^,^ , which 
can be explained if Prsa had a valence of -1-4. Studies 
on the influence of ion-size in Ri_xPrxBa2Cu307_5 show 
that at constant Pr concentration, decreases while the 
antiferromagnetic ordering temperature of Pr ions, Tn, 
increases with the ionic radius of R . It is reasonable 
to assume that the larger the ionic radius of the host rare 
earth R, the easier its substitution by Pr, leading also to 
the higher Tn. Lastly, Zou et al. IJ] found an anoma- 
lously short c lattice parameter in non-superconducting 
PrBa2Cu307_5, which could be due to the smaller size 
of the Pr** ion. 

We have further analyzed the anisotropies shown in 
Fig. |l| and Fig. |^ by computing the power spectral 

density: (-^ioo/oio(p) - Jiio{p)) exp {-ipr} dp . 

Compton profiles, and hence anisotropies, are modu- 
lated with oscillations corresponding to characteristic dis- 
tances over which wavefunctions are coherent in given 
crystallographic directions |^,^. Peaks in the power 
spectral density then indicate these characteristic dis- 
tance scales as shown in Fig. ||. Only the positions and 
relative intensities of the peaks in each spectrum are 
significant since the theoretical anisotropies have been 
scaled. The top panel shows that experiment (solid line 
with dots) and theory (solid line) for YBa2Cu307_5 as 
well as theory for PrBa2Cu307_5 (dashed line) are in 
agreement and that the charge carriers are of a delo- 
calized nature. The periodic structure of the anisotropy 
shown in Fig. |l] and the corresponding peaks in the power 
spectral density in Fig. |3|are directly related to the phases 
of the wavefunction associated with the four O atoms 
surrounding a given Cu atom in the Cu-0 layers [ p9[ . 
In the bottom panel the experiment for PrBa2Cu307_a 
(solid line with dots) shows both a tendency for strong 
suppression of the longer range coherence peaks as in the 
disordered model (sofid fine), as well as a shifting of the 
remaining coherence peak to shorter length scales similar 
to the simple model of insulating, undoped Cu-0 planes 
(dashed line) . This is a clear indication of hole wavefunc- 
tion localization in the Cu-0 layers and of the disorder 



on the Pr and Ba sites. On- going experiments in under- 
doped YBa2Cu307_5 confirm the possibility to monitor 
wavefunction coherence length behaviour near the metal- 
insulator transition using Compton profiles anisotropies. 

With these arguments in mind we can state that super- 
conductivity in PrBa2Cu307_5 appears with decreasing 
PrBa defects and that this is also the reason for the scat- 
ter in the measured Tc values which correspond to vary- 
ing degrees of such disorder, dependent on the method 
and particular conditions of sample fabrication. Mazin 
Q has independently come to a similar conclusion ar- 
guing however that disorder effectively localizes holes in 
the Cu-0 layers in PrBa2Cu307_5 which are of a dif- 
ferent nature with respect to YBa2Cu307_5 due to hy- 
bridization. This possible second mechanism for hole de- 
pletion would also lead to a metal-insulator transition 
in the Cu-0 layers and therefore to similar results for 
the anisotropy. One or both of these mechanisms may 
be active and we point out that disorder-induced dis- 
ruption of charge transfer would provoke insulating be- 
haviour independently of the precise nature of the holes 
which would have been created in the absence of disor- 
der. Further support for this picture is provided by the 
fact that by using a combined strategy of minimizing sub- 
stitution disorder by quenching the melt from high tem- 
peratures after crystal growth and avoiding tetravalent 
PrBa, by using very low oxygen partial pressures during 
growth, we have succeeded in obtaining small flux-grown 
crystallites of PrBa2Cu307_5 with a Tc of 90 K and a 
bulk Meissner effect [^^. The oxygenation used for 
these samples was the same as for 90K superconducting 
YBa2Cu307_5 implying similar, optimal doping in both. 
We expect the electron momentum density in supercon- 
ducting PrBa2Cu307_5 and YBa2Cu307_5 to be similar, 
but at the moment such samples are much too small for 
Compton profile measurements. 

In conclusion we provide a consistent explanation of 
the observed behaviour in PrBa2Cu307_5. Our data 
explicitly shows a strong suppression of the wavefunc- 
tion coherence length in PrBa2Cu307_5 due to disor- 
der in the Ba and Pr planes and the metal-insulator 
transition of the Cu-0 layers. The presence of some 
Pr on the Ba site, probably in a 4+ oxidation state, 
while maintaining the carrier concentration in the Cu-0 
chains, inhibits hole-doping in the Cu-0 planes and/or 
localizes doped carriers. We further propose that su- 
perconductivity appears if this substitution is suppressed 
and our preliminary experiments yielding superconduct- 
ing PrBa2Cu307_5 strongly support this model. 

We acknowledge help from P. Fajardo, D. Vasumathi 
and B. Revaz, and discussions with S. Massidda, S. 
Ishibashi and T. Jarlborg. 
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FIG. 2. Anisotropy (J[ioo/oioi ~ J[iio]) in the electron mo- 
mentum density of PrBa2Cu307_5. The experiment differs 
remarkably from the LDA calculation (Fig. |l|). Two simple 
theoretical models accounting for Praa disorder (solid line) 
and undoped Cu-0 planes (dashed line) are also shown. 



FIG. 3. Power spectral densities of anisotropies. The 
peaks indicate characteristic distance scales over which wave- 
functions are coherent. Top panel: Experiment (solid line 
with dots) and theory (solid line) for YBa2Cu307_i; The- 
ory for PrBa2Cu307_i (dashed line). Bottom panel: In 
PrBa2Cu307_i, experiment (solid line with dots) shows a 
strong suppression of the longer range coherence peaks as in 
the disordered model (solid line) as well as a shifting of the 
remaining coherence peak to shorter length scales similar to 
the model for insulating, undoped Cu-O layers (dashed line). 



FIG. 1. Anisotropy (J[ioo/oioi ~ J[iioi) in the electron mo- 
mentum density of YBa2Cu307_i. Note the good agree- 
ment between the experiment and the LDA calculations 
(YBa2Cu307_i : solid line PrBa2Cu307_5: dashed line). Cal- 
culations are scaled to the experimental variation. 
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